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Moment Analysis and Zipf Law 
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Abstract. The moment analysis method and nuclear Zipf 's law of fragment size distributions are reviewed 
to study nuclear disassembly. In this report, we present a compilation of both theoretical and experimental 
studies on moment analysis and Zipf law performed so far. The relationship of both methods to a possi- 
ble critical behavior or phase transition of nuclear disassembly is discussed. In addition, scaled factorial 
moments and intermittency are reviewed. 
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1 Introduction 

Hot nuclei can be formed in energetic heavy ion colli- 
sions (HIC) and deexcite by different decay modes, such 
as evaporation and multifragmentation. Experimentally, 
multifragment emission was observed to evolve with ex- 
citation energy. The multiplicity, Ni m f, of intermediate 
mass fragment (IMF) rises with the beam energy, reaches 
a maximum, and finally falls to a lower value. The onset of 
multifragmentation may indicate the coexistence of liquid 
and gas phases [1] . Phenomenologically, the mass (charge) 
distribution of IMF distribution can be expressed as a 
power law with parameter T e ft, and a minimum T m ;„ of 
r e ff emerges around the onset point, which suggests that 
a kind of critical behavior may take place. In the frame- 
work of Fisher's droplet model, the mass distribution can 
be described by a power law with a critical exponent of 
r ~ 2.3 when the system is in the vicinity of the critical 
point [2]. 

On the other hand, the caloric curve measurement can 
also provide useful information on the liquid-gas phase 
transition |3|4|5|6|7] , The analysis of other independent 
critical exponents provides additional indications of criti- 
cal behavior of finite nuclear systems |8l9ll0lllll2ll3j . In 
addition, more observables have been proposed to sign the 
liquid-gas phase transition or critical behavior of nuclei 
[14115116117118 19 . Some reviews can be found in this vol- 
ume |20)21|22|23|24] . 

In this report, we shall review the moment analysis 
method and Zipf law of fragment size distribution. The 
phenomenological basis of moment analysis is introduced 
in Sec. 2. Finite size effects are discussed in Sec. 3. Section 4 
gives the application of moment analysis to multifragmen- 
tion and its relation to critical behavior. Scaled factorial 
moments and intermittency are discussed in Sec. 5. In Sec. 



6 Zipf law is introduced for the nuclear fragment distri- 
bution and the corresponding simulations are given; some 
experimental indications of nuclear Zipf law are presented 
in Sec. 7; finally the summary and outlook are given in 
Sec. 8. 



2 Phenomenological Basis of Moment 
Analysis 
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Campi 25 26J and Bauer [27 28] et al. first suggested that 
the methods used in percolation studies may be applied 
to nuclear multifragmentation data. In percolation theory 
the moments of the cluster distribution contain a signature 
of critical behavior [29] . The method of moment analysis 
has been experimentally used to search for evidence of the 
critical behavior in multifragmentation. The definition of 
the k moments of the cluster size distribution for each 
event is 

M k = £ A k n A , (1) 

where A is the fragment mass, and tia is the number of 
charged fragments whose charge is Z and mass is A. The 
sum runs over all masses A in the event including neutrons 
except the heaviest fragment {A max ). This quantity was 
taken as a basic tool in extracting critical exponents in 
Au + C data [9] . It has been argued that there should be 
an enhancement in the critical region of the moment M/-, 
for k > t — 1 , with a critical exponent r > 2 [25126] . 

In experimental analyses, events are sorted by differ- 
ent conditions. In this case, so-called conditional moments 
are used to describe the fragment distribution. Usually 
the mean value of Mk{m) for events with given control 
parameters, e.g. the moment Mfc for events with a given 
multiplicity to, or total bound charge number Zb oun d, or 
excitation energy E* , is called conditional moment. 
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More insight in the shape of the fragment size distribu- 
tion is obtained by looking at a combination of moments 
Mk- For example, the quantity 



72 



M 2 M 
Ml 



1, 



(2) 



has been used, where Mi and M 2 are the first and sec- 
ond moments of the mass distribution and Mo is the total 
multiplicity including neutrons, a 2 is the variance of the 
fragment distribution and (s) = Mi /Mo represents the 
mean fragment size. 72 takes the value 72 = 2 for a pure 
exponential distribution N(s) ~ exp(-as) regardless of 
the value of a, but 72 3> 2 for a power-law distribution 
N(s) ~ s~ r when r > 2. In the percolation model, the po- 
sition of the maximum 72 value defines the critical point, 
where the fluctuations in the fragment size distribution 
are the largest. In principle a genuine critical behavior re- 
quires the peak value of 72 to be larger than 2 |25|26j . 
However, due to finite size effects, this is not always true 
when the system size decreases, as we will see in the fol- 
lowing sections. 

Campi also suggested to use the single event (j) mo- 
ment, i.e. 

(3) 



K = 



A^A mc 



to investigate the shape of fragment size distribution. Also 
normalized moments 1251 



Sjp=Mg>/M? 



r(j) /i\Ai) 



(4) 



can be defined. It was suggested to use the event-by-event 
scatter-plots of the natural log of the size (A max ) or charge 
number (Z max ) of the largest cluster, lnA max or h\Z max 
versus the natural log of the second moment, lnM 2 , or the 
normalized moment lnS^ to search for the largest fluctu- 
ation point. Some examples will be given in the following 
sections. 

In the percolation model, the cluster size distribution 
for infinite systems near a critical point can be expressed 

by 



n{s) ~ s- T f{es°). 



(5) 



where s is the size of finite clusters, r and a two critical 
exponents and e a variable that characterizes the state of 
the system. In thermal phase transitions e = T — T c is the 
distance to the critical temperature T c . In percolation e = 
p c —p is the distance to the critical fraction of active bonds 
or occupied sites p c . The scaling function /(es CT ) satisfies 
/(0) = 1, decaying rapidly (exponentially) for large values 
of |e|. In addition, theory predicts that when e < one 
infinite cluster (liquid or gel) is present in the system while 
no such cluster exists when e > (only droplets or Tu- 
rners) . In finite systems a similar behavior is observed, 
especially when the largest cluster is counted separately. 

The moment analysis method is useful to obtain some 
information about the possible occurrence of a critical be- 
havior. In general, critical exponents can be defined ac- 
cording to the standard procedure followed in condensed 



matter physics [3U]. For example, 
M fc (e) = J2 = A k n A (e) ~ 



S(e-O) (6) 



where r and a are the critical exponents. For the percola- 
tion phase transition and the critical point in the Fisher 
droplet model, the exponent r satisfies 2 < r < 3 and thus 
the second and high moments diverge at the critical point. 
In contrast the lower moments Mo and Mi, which corre- 
spond to the number of fragments and the total mass, do 
not diverge. 

Based upon the scaling relation Eq. (JSJ), there exists 
the following relationship between critical exponents and 
moments: 



M ~ 
Mi 

M 2 n 



2-a 



(7) 



where [3 and 7 are two other critical exponents. Some re- 
lationships among critical exponents exist (hyperscaling 
relations), for instance 



2/3 + 7 



1 



(8) 



In finite systems transitions are smooth, but it is still 
possible to determine some critical exponents, as we will 
discuss in the next section. By analogy with the infinite 
system behavior, one says that these moments exhibit a 
critical behavior also for finite systems. In particular in the 
Fisher model, the thermal critical point is also a critical 
point for moments of the fragment size distribution. 

In order to illustrate the application of moment analy- 
sis, we show the EOS data and NIMROD data as examples 
in Sec. 4. 



3 Finite Size Effects 

Since the nucleus is a finite size system, the macroscopic 
thermal limit cannot be applied. Therefore finite size ef- 
fects on phase transition behavior should be checked. In 
this section, we give some examples to illustrate this prob- 
lem. 

A percolation on a cubic lattice of linear size L contain- 
ing L 3 sites, for L — 4 to 10, where all sites are occupied 
and bonds are assumed to exist between neighbouring sites 
with bond probability p, has been considered [31]. Sites 
that are connected together by such bonds are said to be- 
long to the same cluster. It is well known that in such a 
model there exists a critical (or threshold) probability p c 
such that for p > p c there is a large cluster that percolates 
throughout the lattice from end to end whereas for p < p c 
no such cluster exists and all the sites belong to small 
clusters (including isolated sites, i.e. singlet or clusters of 
size 1). As L — > 00 the transition becomes sharper and 
p c approaches a limiting value which for bond percolation 
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on a cubic lattice is p c = 0.249 [29]. For finite systems the 
threshold percolation probability is not so sharply defined. 

In order to quantitatively illustrate finite size effects 
on critical behavior, the average of normalized second mo- 
ment 52 (£2) over all events belonging to the same value 
of \n(A max ) was calculated [31]. The results obtained by 
such averaging are presented by the dots shown in Fig. [1] 
for various cubic lattices with linear dimension L = 4 - 
10 sites [31]. The location of the maximum value of S2 is 
now defined as corresponding to the location of the critical 
point, which is a standard way of determining the perco- 
lation threshold [32]. The slope of the lower branches of 
the curves in Fig. Q] can also be calculated. This slope is 
expected to be 1 + which for percolation in three di- 
mensions is equal to 1.23. For comparison the slopes of the 
straight lines by a lest-squares fit to the lower branches 
of the L — 4 to 10 curves are found, in ascending or- 
der of L, to have the values 1.582±0.036, 1.503±0.029, 
1.375±0.017, 1.355±0.021, 1.260±0.007, 1.258±0.014 and 
1.242±0.015 [31 . This indicates that these slopes rapidly 
approach the value expected in the thermodynamic limit. 
In calculating these slopes one has excluded the points 
near the bottom of the branch in the region where the 
curves in Fig. [1] deviate noticeably from a straight line. 
These points correspond to events that are far from the 
critical region. 

Similarly to the analysis for the correlation of S2 and 
\n(A max ), finite size effects have been also investigated for 
Mi by Campi [26]. This is shown in Fig. [2] where M2(n) 
is plotted for various system sizes (50 3 , 9 3 , 5 3 and 3 3 ) in a 
percolation model. We see clearly the critical behavior for 
the largest system, namely a well defined peak, and how 
this peak is smoothed when decreasing the size |26j . 

4 Application of the Moment analysis method 
4.1 EOS data 

4.1.1 Experimental description 

The reverse kinematic EOS experiment was performed 
with 1 A GeV 197 Au, 139 La, and 84 Kr beams on carbon 
targets. The experiment was done with the EOS Time 
Project Chamber (TPC) and multiple sampling ionization 
chamber (MUSIC II). The excellent charge resolution of 
this detector permitted identification of all detected frag- 
ments. The fully reconstructed multifragmentation events 
for which the total charge of the system was taken as 
79 < Z < 83, 54 < Z < 60, 3 3 < Z < 39 for Au, La, and 
Kr, respectively [33 34 35 36] were analyzed. The remnant 
refers to the equilibrated nucleus formed after the emis- 
sion of prompt particles. The charge and mass of the rem- 
nant were obtained by removing for each event the total 
charge of the prompt particles. The excitation energy of 
the remnant E* was based on an energy balance between 
the excited remnant and the final stage of the fragments 
for each event [37] • The thermal excitation energy E^ h of 
the remnant was obtained as the difference between E* 
and E x which is a nonthermal component, namely an ex- 
pansion energy [33 34 35 36138] . 
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Fig. 1. The logarithm of the largest fragment size A max as a 
function of the logarithm of the corresponding average normal- 
ized second moment S2 for bond percolation on simple cubic 
lattices of linear size ranging from L = 4 - 10 sites. The dots 
represent the actual calculation results and the curves drawn 
are just to guide the eye. The number next to each curve gives 
the value of the linear size L. Note that the ln52 scale given 
corresponds to the L = 10 curve. The other curves are succes- 
sively shifted to the left with respect to each other by a distance 
of 0.25. The dashed curve and the dotted-dashed straight line 
are explained in the text. Figure is taken from Ref. [31] . 

4.1.2 Determination of Critical Point and Exponent in 
Terms of Moment Analysis 

The determination of the critical point and the associ- 
ated exponents in the multifragmentation of gold nuclei 
was first attempted by the EOS collaboration [9i[. In their 
early publication [9], they use the multiplicity m, as a con- 
trol variable for the collision violence and assume that m 
is a linear measure of the distance from the critical point. 
Then the critical exponents /?, 7 and r, can be determined 
according to Eqs. (7,8) above. They find that these expo- 
nents are close to the nominal liquid-gas universality class 
values. However, this method is very delicate. In partic- 
ular, due to the small size of the system, an important 
rounding of the transition is expected which may distort 
considerably the determined critical exponents. For a re- 
view of this debate, see the arguments between Bauer [39] 
and Gilkes [40] . 

A different analysis was also proposed by the EOS col- 
laboration [41] . In this work, thermal excitation energy 
has been taken as a control variable, which is believed to 
be more suitable to characterize the collision violence. 

The 72 analysis is shown in Fig. [3] for all three systems. 
The position of the maximum 72 value defines the critical 
excitation energy E*, which corresponds to the largest 
fluctuation point in the fragment size distribution. The 
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Fig. 2. The conditional moments Ma(n) for percolation in a 
cubic lattice of linear size L = 3, 5, 9 and 50 (the corresponding 
cubic lattice is L z which is shown as the number in the insert). 
The figure is taken from Ref. [26]. 



peak in 72 is well denned for La and Au. For Kr, the peak 
is very broad and the value 72 is less than 2. 

Fig. [3] also shows a 72 calculation using the statisti- 
cal multifragmentation model (SMM). The fission contri- 
bution to 72 has been removed both from the data and 
SMM. In the case of Au, the 72 value remains above two 
for most of the excitation energy range both in data and 
SMM. The E* h width over which 72 > 2 is smaller for La 
and disappears for Kr. The decrease in 72 with decreas- 
ing system size is also seen in 3D percolation studies and 
these differences have been attributed to finite size effects 
[41142143] . 

The exponent r can be obtained if the second moment 
M2 and the third moment M3 of the fragment mass dis- 
tributions are known. A plot of ln(Ma) vs ln(M2) should 
give a straight line with a slope given by 



Aln(M 3 ) 
Aln{M 2 ) 



t-3 



(9) 



Fig. Q] shows a scatter-plot of ln(M3) vs ln(M2) for the 
three systems constructed with data above the critical ex- 
citation energy E* (see Fig. [3]) and with SMM simulations. 
A linear fit to \sx{M^) vs ln(M2) gives the value of r. The 
fitted t values are 2.16 ± 0.08, 2.10 ± 0.06 and 1.88 ± 0.08, 
respectively. The former two are very close to the critical 
exponents r ~ 2.3 of the liquid-gas universal class. 

The exponent (3 can be obtained for the multifragmen- 
tation data by the relation 
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Fig. 3. 72 as a function of E^ h for all three systems of 1 A GeV 
Au, La, and Kr collisions with C target and SMM calculations. 
Figure is taken from Ref. [33] . 
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Fig. 4. ln(Ms) vs ln(M2) for Au, La, and Kr above the critical 
energy. Figure is taken from Ref. [33] . 
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Fig. 5. \n(A max ) vs \n\El h — Eq\ for Au, La, and Kr below the 
critical energy for exponent (3 determination. Figure is taken 
from Ref. 1331. 



Fig. 6. Scatter-plots of \n(A max ) vs ln(M2) from the data 
for Au, La, and Kr. Left panel: EOS data; Right panel: SMM 
simulation. The figure is taken from Ref. [33]. 



where e = p — p c and e > 0. In the multifragmentation 
case of this work p and p c have been replaced by E^ h and 
E*. In an infinite system, the finite cluster exists only on 
the liquid side of p c . In a finite system a largest cluster is 
present on both sides of the critical point, but the above 
equation holds only on the liquid side. Fig. [5] shows a plot 
of \n(A max ) vs \n\Et h — E*\ for Au, La, and Kr. The values 
of (3 extracted for Au and La are 0.32±0.02 and 0.34±0.02, 
respectively, which are close to the value of 0.33 predicted 
for a liquid-gas phase transition. On the other hand, the 
value of j3 = 0.53 ± 0.05 for Kr is much higher than that 
of Au and La. 

As shown in Sec. 2, Campi also suggested that the cor- 
relation between the size of the biggest fragment A max and 
the moments in each event, i.e. the scatter-plot, can mea- 
sure the critical behavior in nuclei. Fig. O depicts a scatter- 
plot with logarithmic scale for Au, La, and Kr of EOS 
data. The two branches corresponding to the sub-critical 
(upper branch) and overcritical (lower branch) events are 
clearly seen for Au and La. The scatter-plot is very broad 
for Kr and fills most of the available space. The sub- and 
over-critical branches seem to overlap and are not well sep- 
arated. Studies on SMM show a similar behavior. If one 
knows the location of the critical point from some other 
methods, then the scatter-plot can be used to calculate the 
ratio of critical exponents from the slope of the sub- 
critical branch. In EOS data, the position of the largest 72 
was used to define the critical point, which corresponds to 
the largest fluctuation of the fragment distribution. In this 
context, (3/j values for Au, La and Kr can be extracted 
from the linear fit to the upper branch; they are 0.22 ± 



0.03, 0.25 ± 0.01 and 0.50 ± 0.01, respectively. values 
of Au and La are close to the value 0.26 expected for the 
liquid-gas universality class. 

To summarize the critical exponent analysis of the 
EOS data, the experimental results in conjunction with 
SMM provide some indications on the order of the phase 
transition in Au, La and Kr. The values of the critical 
exponents r, /3, and 7, which are close to the values of a 
liquid-gas system, along with nearly zero latent heat (this 
subject is beyond the discussion topics in this review, but 
the interested reader is reported to Refs. [33134] ) have been 
interpreted by the authors as a continuous phase transi- 
tion in Au and La. However, the analysis of Kr leads to 
very different critical exponents. A recent analysis based 
on the shape of SMM microcanonical caloric curve indi- 
cates a first order phase transition for the multifragmen- 
tation of Kr [33134] . 

4.2 NIMROD data 

4.2.1 Experimental set-up and Analysis Details 

Using the TAMU NIMROD (Neutron Ion Multidetector 
for Reaction Oriented Dynamics) and beams from the 
TAMU K500 super-conducting cyclotron, we have probed 
the properties of excited projectile- like fragments produced 
in the reactions of 47 MeV/nucleon 40 Ar + 27 Al, 48 Ti and 
58 Ni. The charged particle detector array of NIMROD, 
which is set inside a neutron ball, includes 166 individual 
Csl detectors arranged in 12 rings in polar angles from 
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Fig. 7. Campi plots for nine intervals of excitation energy for 
the QP formed in 40 Ar + 58 Ni. Figure is taken from Ref. [14]. 



~ 3° to ~ 170°. The detailed description for the experi- 
ment can be found in |14j . The correlation of the charged 
particle multiplicity (M cp ) and the neutron multiplicity 
(M n ) was used to sort the event violence. After the recon- 
struction of the quasi-projectile (QP) particle source, the 
excitation energy was deduced event-by-event using the 
energy balance equation [37] , 



4.2.2 Critical Point Determination via Moment Analysis 

In Fig. [7J we present Campi scatter-plots for the nine se- 
lected excitation energy bins. In the low excitation en- 
ergy bins of E*/A < 3.7 MeV/u, the upper (liquid phase) 
branch is strongly dominant while at E* j A > 7.5 MeV/u, 
the lower Z max (gas phase) branch is strongly dominant. 
In the region of intermediate E*/A of 4.6- 6.5 MeV/u, the 
transition from the liquid dominated branch to the va- 
por branch occurs, indicating that the region of maximal 
fluctuations is to be found in that range. 

The excitation energy dependence of the average val- 
ues of 72 obtained in an event-by-event analysis of our 
data are shown in Fig. [8] 72 reaches its maximum in the 
5-6 MeV excitation energy range. In contrast to observa- 
tions for heavier systems of Au and La [33141] , there is no 
well defined peak in 72 for our very light system and 72 
is relatively constant at higher excitation energies. This is 
similar to the case of Kr of EOS data. We note also that 
the peak value of 72 is lower than 2 which is the expected 
smallest value for critical behavior in large systems. How- 
ever, 3D percolation studies indicate that finite size effects 
can lead to a decrease of 72 with system size [42143] . For 
a percolation system with 64 sites, peaks in 72 under two 
are observed. Therefore, the criterion 72 > 2 alone is not 
sufficient to discriminate whether or not the critical point 
is reached. 
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Fig. 8. 72 of the QP systems formed in Ar + Al (open cir- 
cles), Ti (open triangles) and Ni (solid squares) as a function 
of excitation energy. 
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Fig. 9. Charge distribution of QP in different E* /A window 
for the reaction 40 Ar + 58 Ni. Lines represent fits. Figure is 
taken from Ref. I14I. 



In the Fisher droplet model, the critical exponent t 
can be deduced from the cluster distribution near the crit- 
ical point. To quantitatively pin down the possible phase 
transition point, we use a power law fit to the QP charge 
distribution in the range of Z — 2 - 7 (Fig. [5]) to extract 
the effective Fisher-law parameter r e ff by 



dN/dZ - Z- T '" . 



(11) 



Fig- flOT a) shows the effective Fisher- law parameter r e f f 
as a function of excitation energy. A minimum with r e // 
~ 2.3 is seen to occur in the E* /A range of 5 to 6 MeV/u 
[44] . This value is close to the critical exponent of the 
liquid-gas phase transition universality class [2]. 

Assuming that the heaviest cluster in each event rep- 
resents the liquid phase, we have attempted to isolate the 
gas phase by event-by-event removal of the heaviest clus- 
ter from the charge distributions. We find that the resul- 
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Fig. 10. The effective Fisher-law parameter (r e ff) (a), the 
effective exponential law parameter (A e //) (b), (S2) (c), NVZ 
fluctuation (d), the mean charge number of the second largest 
fragment {Z2max) (e), the Zipf-law parameter A (f). See details 
in text. Figure is taken from Ref. |44] , 



tant distributions are better described with an exponential 
form exp~ XcffZ . The fitting parameter X e ff was derived 
and is plotted against excitation energy in Fig. [TOT b). A 
minimum is seen in the same region where r e ff shows a 
minimum. To further explore this region we have inves- 
tigated other proposed observables commonly related to 
fluctuations and critical behavior. Fig. [TOf c) shows the 
mean normalized second moment, {S2) as a function of 
excitation energy. A peak is seen around 5.6 MeV/u, it 
indicates that the fluctuation of the fragment distribu- 
tion is the largest in this excitation energy region. Simi- 
larly, the normalized variance in Z max /ZQp distribution 

(i.e. NVZ = {z" 1 "* /Zqp) ) S3 shows a maximum in the 
same excitation energy region (Fig. llOf d)). which illus- 
trates the maximal fluctuation for the largest fragment 
is reached around E*/A = 5.6 MeV. The second largest 
fragment shows a behavior similar to the one of the largest 
fragment. Fig. [T07 e) shows a broad peak of (Z2 ma x) - the 
average atomic number of the second largest fragment - 
also occurring in the same excitation energy range around 
5.6 MeV/u. 

More variables have been collected to support the de- 
termination of the critical point around 5.6 MeV/u of ex- 
citation energy for our system [14] , such as Z\-scaling [46] 
or energy fluctuations 23J. In addition, the measurement 
of the caloric curve [14] gives a temperature, T c ~ 8.3 
MeV around E*/A = 5.6 MeV. The value of the critical 
temperature is needed for the determination of the critical 
exponents, as explained in the following subsection. 
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Fig. 11. The correlation between In (S3) vs ln(5*2) and a linear 
fit. 



4.2.3 Determination of Critical Exponents Based on 
Moment Analysis 

In terms of the scaling theory, t can also be deduced from 
Eq. ©. Since the value of T c = 8.3 MeV has been deter- 
mined from our caloric curve measurements [14] . we can 
explore the correlation of S2 and S3 in two ranges of ex- 
citation energy (see Figure [TTj) . The moments were calcu- 
lated by excluding the species with Z max for the "liquid'' 
phase but including it in the "vapor" phase. The slopes 
were determined from linear fits to the "vapor" and "liq- 
uid" regions respectively and then averaged. In this way, 
we obtained a value of r = 2.13 ± 0.1. 

Other critical exponents can also be related to other 
moments of the cluster distribution, M^, Using our caloric 
curve measurements [14] , we can use temperature as a con- 
trol parameter for such determinations. Then the critical 
exponent j3 can be extracted from the relation 



oc (1 



T> ' 



(12) 



and the critical exponent 7 can be extracted from the 
second moment via 



M 2 oc|i-^|-r 



(13) 



In both equations, |1 — is the parameter which mea- 
sures the distance from the critical point. 

The upper panel of Fig. [T2] explores the dependence 
of Z max on (1 — ;yr). A dramatic change of Z max around 
the critical temperature T c is observed. Lattice-gas model 
(LGM) calculations also predict that the slope of Z max 
vs T will change at the liquid-gas phase transition [47] . 
Using the liquid side points, we can deduce the critical 
exponent (3 by ln(Z max ) vs ln\l — T/T c \. Fig.fLWa) shows 
the extraction of (3 using Eq. (fL2|l . An excellent fit was 
obtained in the region away from the critical point, which 
indicates a critical exponent (3 = 0.33 ± 0.01. Near the 
critical point, finite size effects become stronger so that 
the scaling law is violated. The extracted value of (3 is that 
expected for a liquid-gas transition (See Table 1) [29] . 
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Table 1. Comparison of the Critical Exponents 
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Fig. 12. The extraction of the critical exponent j3 (a) and 7 
(b). See texts for details. 



To extract the critical exponent 7, we take M2 on the 
liquid side without Z max . Fig. [T^b) shows ln(M2) as a 
function of ln( 1 1 — ^- | ) . We center our fit to Eq. (TP3"]) about 

the center of the range of (1 — T/T c ) which leads to the 
linear fit and extraction of (3 as represented in Figure [P2l 
We obtain a critical exponent 7 = 1.15 ± 0.06. This value 
of 7 is also close to the value expected for the liquid-gas 
universality class (see Table 1). It is seen that the selected 
region has a good power law dependence. 

Since we have the critical exponent (3 and 7, we can 
use the scaling relation 



1 



P + l 



(14) 



to derive the critical exponent a. In such way, we get the 
g = 0.68 ± 0.04, which is also very close to the expected 
critical exponent of a liquid-gas system. 

To summarize the critical exponents extracted from 
NIMROD data, we present the results in Table 1 as well 
as the values expected for the 3D percolation and liquid- 
gas universality classes. It is apparent that our values for 
this light system with A~36 are closer to the values of the 
liquid-gas phase transition universality class rather than 
to the 3D percolation class. 



Exponents 


3D Percolation 


Liquid-Gas 


NIMROD 


T 


2.18 


2.21 


2.13L0.10 


P 


0.41 


0.33 


0.33L0.01 


7 


1.8 


1.23 


1.15L0.06 


a 


0.45 


0.64 


0.68±0.04 



5 Scaled factorial moments and intermittency 

Intermittency is related to the existence of large non- 
statistical fluctuations and is a signal of self-similarity of 
the fluctuation distribution at all scales. This signal can 
be deduced from the scaled factorial moments 051, 



F k (6) 



c/5 



(ni(rii - l)(n, - 2)...{rn -k + l)} 



i=l 



c/6 



(rn) k 



(15) 

where X max is an upper characteristic value of the system 
(i.e. total mass or charge, maximum transverse energy or 
momentum, etc.) and k is the order of the moment. The 
total interval 0-X max (1-A max , Z max in the case of mass 
or charge distributions) is divided into X max /8 bins of 
the size 5, rii is the number of particles in the ith bin 
for an event, and the ensemble average () is performed 
over all events. The concept of intermittency was origi- 
nally developed in the field of fluid dynamics to study the 
fluctuations occurring in turbulent flows [49150] . Its pres- 
ence in the velocity and temperature distributions is estab- 
lished by the existence of large non-statistical fluctuations 
which exhibit scale invariance. Intermittency in physical 
systems is studied by examining the scaling properties of 
the moments of the distributions of relevant variables over 
a range of scales [51]. The concept of intermittency was 
first introduced for the study of dynamical fluctuations in 
the density distribution of particles produced in high en- 
ergy collisions by Bialas and Peschanski [48] . It soon led to 
the discovery of a characteristic power law dependence of 
the factorial moments, F kl of an order k on the resolution 
scale, 5: F k oc (1/S)^ k \ The specific properties of the in- 
termittency exponent, /(fc), can be associated either with 
a random production process [48l52j or with a second- 
order phase transition [52153154] depending on the values 
obtained. Thus an analysis of the factorial moments may 
provide important information on the dynamical proper- 
ties of the system. Ploszajczak and Tucholski were the first 
to suggest searching for intermittency patterns in the mass 
and charge distributions of the fragments produced in en- 
ergetic collisions [55]. Since then many studies show that 
an intermittency pattern of fluctuations in the fragmenta- 
tion charge distributions has been observed in many data 
and models. Much effort has been devoted to find the re- 
lation between fragmentation, a possible critical behavior, 
and intermittency [45|56|57|58|59|60] . 

Intermittency is defined by the relation 



F k (5') = F k (aS)=a-f^F k (S), 



(16) 
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between factorial moments Fk(S') and Fk(S s ) obtained for 
two different binning parameters S and 5' = a<5. Intermit- 
tency implies a linear relationship in the double logarith- 
mic plot of InFk versus —InS. 

The fractal intermittency exponent, f(k), is related to 
the factorial dimension dk by 



/(*) 



k - 1 



> 0. 



(17) 



Different processes seem to give a different behavior of 
these anomalous fractal dimension dk'- (1) dk = constant 
corresponds to a monofractal, second order phase transi- 
tion in the Ising model and in the Feynman- Wilson fluid 
[53.55- ^ nas been also demonstrated that in the case of a 
second order phase transition in the Ginzburg-Landau de- 
scription one gets dk = d%(k — l)' t_1 with \i = 1.304 [53]. 
(2) dk oc k correspond to multifractal, cascading processes 
[48] , Therefore, a study of the anomalous fractal dimen- 
sions can give useful information about the evolution of 
the system. 

Several models have been introduced to study the in- 
termittency signal. One of the simplest models, widely 
used in the analysis of experimental data and which gives 
intermittency, is the percolation model. Percolation mod- 
els predict a phase transition corrected for finite size ef- 
fects and produce, at the critical point for this phase tran- 
sition, a mass distribution following a power law and obey- 
ing scaling properties. 

An intermittency analysis has been performed on many 
heavy ion collision data as well as emulsion data. Here we 
give an example of the multifragmentation data of Au + 
Au collisions at 35 MeV/u which was performed at NSCL 
by the Multics-Miniball Collaboration [61 . A power-law 
charge distribution, A~ T with r ~ 2.2 and an intermit- 
tency signal has been observed for the events selected in 
the region of the Campi scatter-plot where "critical" be- 
havior is expected. As shown in Fig. [T3J three cuts have 
been tested. The upper branch is mostly related to the liq- 
uid branch and the lower branch to the gas branch, while 
the central cut (2) is expected to belong to a region where 
critical behavior takes place. Actually the resultant charge 
distribution of cut (2) shows a power-law distribution with 
r ~ 2.2 which is close to the droplet model prediction if 
the liquid-gas critical point is explored. The scaled facto- 
rial moments are shown in Fig. 14 for the different cuts 
of Figure 13. For cut 3, the logarithm of the scaled fac- 
torial moment is always negative and almost independent 
of -ln<5; there is no intermittency signal. The situation is 
different for cut 2 (the central part). The logarithm of the 
scaled factorial moments is positive and almost linearly 
increasing as a function of -ln<5, and an intermittency has 
been observed. Cut 1 gives a zero slope, no intermittency 
signal again. 

It has been argued that the interpretation of this ex- 
perimentally observed intermittency signal may, however, 
be problematic due to an ensemble average effect [55] . 
Since cut 2 involves a large range of impact parameters, 
the observed intermittency signal could be an artifact of 
ensemble averaging, and can not be seen as a definite ev- 
idence of large fluctuation driven by a critical behavior. 
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Fig. 13. Experimental Campi-scatter-plots from Ref. [61j . 
Three cnts are employed to selected the upper branch (1), the 
lower branch (3), and the central region (2). 
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Fig. 14. Experimental results from Ref. [61]. Scaled factorial 
moments \n(Fk) vs -ln(5 3 ) for the three cuts made on Fig. [13] 
left part cut 1, central part cut 2, and right part cut 3. Solid 
circles represent the SFM of order k = 2, open circles k = 3, 
open squares k = 4, and open triangles k = 5. Figure is taken 
from Ref. [61] . 



Actually, several criticisms have been raised about the 
role of the intermittency signal in nuclear fragmentation. 
For instance, Elattari et al. showed that an intermittency 
signal can be obtained even for a simple fragmentation 
generator model by the random population of mass bins 
with a power law distribution in which the only non- 
statistical source of fluctuations is the mass conservation 
law [57]- It has also been shown that the intermittency 
signal is washed out when events of fixed total multiplic- 
ity are selected [45|60j or when the size of the system 
tends to infinity in the percolation model in which the 
fluctuations are of nontrivial origin [60]. Moreover, the 
intermittency signal is not observed in the narrow exci- 
tation energy region where the phase transition occurs in 
the framework of the well-known Copenhagen statistical 
multi-fragmentation model [58] or in the data of 35-110 
MeV/nucleon 36 Ar+ 197 Au when the effects of impact pa- 
rameter averaging are reduced by some appropriate cuts 
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Fig. 15. Left panels: the multiplicity distribution (upper 
panel), the mass distribution (middle panel), the scaled facto- 
rial moments (bottom panel) with the multiplicity restriction 
for 129 Xe in the lattice gas model calculation. Right panel: 
Same as the left panel but for the events mixed with T — 5.5 
MeV and T = 7 MeV. Figure is taken from Ref. [62]. 



|56j . However, it is important to notice that there is no 
reason to expect intermittency if the phase transition is 
first order. 

As an example, we check the intermittency behavior 
[62] in the Lattice Gas model for the disassembly of the 
system 129 Xe at 0.38/5o in the framework of LGM (for 
the details of the model description, please see the follow- 
ing section). At a temperature T = 5.5 MeV, the mass 
distribution shows a power-law distribution with an effec- 
tive power-law parameter r = 2.43. In a previous work 
with the same model, it was shown that the liquid-gas 
phase transition occurs near 5.5 MeV for this system in 
the LGM |63|64| . The InF^ shows slight negative values 
with slightly positive slopes versus —Ind. However, this 
kind of the positive slopes with a moment less than unity 
may be of trivial origin and does not demonstrate the ap- 
pearance of intermittency which is characteristic of sys- 
tems exhibiting larger than Poisson fluctuations (i.e. the 
moment should be larger than unity). In order to check 
the event mixture effect on the scaled factorial moment, 
we mixed all the events at T = 4 MeV and T = 7 MeV 
and also used the multiplicity cuts (29 < M < 101) and 
(M < 29 or M > 101) to see if an intermittency behav- 
ior can be found in such mixed events. Figure [15] shows 
these results. Even though all the InFk values are positive, 
they are flat, i.e. there is no intermittency signal. In these 
cases, the fluctuation is large enough but the mass distri- 
bution shows no power-law distribution. Hence, intermit- 
tency is absent. However, intermittency emerges when the 



moments were calculated from the mixed events of T = 5.5 
MeV and T = 7 MeV (Fig. [15]). In this case, the mass dis- 
tribution shows a quite good power-law distribution and 
fluctuations are also large enough to induce intermittency. 

From the above discussions, the apparent signals of in- 
termittency which emerge in many experimental data are 
not easy to understand since many experimental condi- 
tions bring some complexities to the pure signal of inter- 
mittency, such as event mixing. More precise experimental 
measurements in the future are needed to probe the inter- 
mittency signal, which then may be taken as a signal of 
true critical behavior. 

6 Phenomenological Basis of Nuclear Zipf 
Law and Model Simulation 

In the above sections, we have focussed on the moment 
analysis, namely the behavior of the moments of the frag- 
ment size distribution, or of the scaled factorial moments. 
Both are related to the fluctuations of some physical ob- 
servables. In this section, we would like to emphasize the 
topological structure of the fragment size distribution, i.e. 
how the fragments distribute from the largest to the small- 
est in nuclear fragmentation. To this end, we introduce the 
Zipf-type plot, i.e. rank-ordering plot, in the fragment size 
distribution as well as Zipf 's law which will be illustrated 
in the following [64|65j . 

The original Zipf's law [66] has been used for the di- 
agnosis of nuclear liquid-gas phase transition and as such 
we have called it the nuclear Zipf's law. Zipf's law has 
been known as a statistical phenomenon concerning the 
relation between English words and their frequency in lit- 
erature in the field of linguistics [53]. The law states that, 
when we list the words in the order of decreasing popu- 
lation, the frequency of a word is inversely proportional 
to its rank [66j . This relation was found not only in lin- 
guistics but also in other fields of sciences. For instance, 
the law appeared in distributions of populations in cities, 
distributions of income of corporations, distributions of 
areas of lakes and cluster-size distribution in percolation 
processes |67l68j . The details for the proposal of nuclear 
Zipf's law can been found in Ref. [64165] . In this report, we 
firstly define the nuclear Zipf plot for the fragment mass 
(charge) distribution and nuclear Zipf's law in the simu- 
lation with help of the lattice gas model. Then we show 
some experimental evidences for the nuclear Zipf law as 
well as some remarks. 

The tools we will use here are the isospin dependent 
lattice gas model (LGM) and molecular dynamical model 
(MD). The lattice gas model was developed to describe the 
liquid-gas phase transition for atomic systems by Lee and 
Yang [55] • The same model has already been applied to nu- 
clear physics for isospin symmetrical systems in the grand- 
canonical ensemble |70j with a sampling of the canonical 
ensemble |63I71I72I73I74I75I76| , and also for isospin asym- 
metrical nuclear matter in the mean field approximation 
[77j . In addition, a classical molecular dynamical model is 
used to compare its results with the results of the lattice 
gas model. 
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Fig. 16. Effective power-law parameter, r, second moment of 
the cluster distribution, S2, and multiplicity of intermediate 
mass fragments, N im f as a function of temperature for the 
disassembly of 129 Xe at p f ~ 0.38p in I-LGM. The arrow 
represents the estimated temperature of the phase transition. 
Figure is taken from 65 . 
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In the lattice gas model, A (= TV + Z) nucleons with 
an occupation number s which is defined s = 1 (-1) for a 
proton (neutron) or s = for a vacancy, are placed on the 
L sites of the lattice. Nucleons in the nearest neighbor- 
ing sites interact with an energy e SiSj . The hamiltonian 

is written as E = J2t=i^ ~ J2i< 3 es iSj SiSj. A three- 
dimension cubic lattice with L sites is used. The freeze- 
out density of disassembling system is assumed to be pf 
= j- po , where po is the normal nuclear density. The disas- 
sembly of the system is to be calculated at pf, beyond 
which nucleons are too far apart to interact. Nucleons 
are put into lattice by Monte Carlo Metropolis sampling. 
Once the nucleons have been placed we also ascribe to 
each of them a momentum by Monte Carlo samplings 
of a Maxwell-Boltzmann distribution. Once this is done 
the LGM immediately gives the cluster distribution us- 
ing the rule that two nucleons are part of the same clus- 
ter if P,?/2yU — e SiSj SiSj < 0. This method is similar to 
the Coniglio-Klein prescription [7R in condensed matter 
physics and was shown to be valid in LGM [71172174175] . In 
addition, to calculate clusters using MD we propagate the 
particles from the initial configuration for a long time un- 
der the influence of the chosen force. The form of the force 
is chosen to compare with the results of LGM. The sys- 
tem evolves with the potential. At asymptotic times the 
clusters are easily recognized. Obscrvables based on the 
cluster distribution in both models can now be compared. 
In the case of proton-proton interactions, the Coulomb 
interaction can also be added separately and it can be 
compared with the case without Coulomb effects. 

In order to check the phase transition behavior in the 
I-LGM, we will first show the calculations of some physical 
observables in Fig. [TH namely the effective power-law pa- 
rameter, r, the second moment of the cluster distribution, 



Fig. 17. Average charge Z n with rank n as a function of n for 
129 Xe pf ~ 0.38po in I-LGM. The histograms are the calcula- 
tion results and the straight lines are their fits with 
Figure is taken from [65] . 



S2 [60] . and the multiplicity of intermediate mass frag- 
ments, Nimf, for the disassembly of 129 Xe at the freeze-out 
density pt ~ 0.38,00 . These observables have been success- 
fully employed in previous works to probe the liquid- gas 
phase transition, as shown in Rcfs. 63 65 75]- The val- 
ley of r, the peaks of A^ m / and S2 are located around 
T ~ 5.5 MeV which is the signature of the onset of the 
phase transition. Because of the exact mapping between 
the LGM and the Ising model, we know that at this point 
the transition is first order. 

Now we present the results for testing Zipf 's law in the 
charge distribution of clusters. The law states that the re- 
lation between the sizes and their ranks is described by 
Z n = c/n (n=l, 2, 3, ...), where c is a constant and Z n (or 
A n ) is the average charge (or mass) of rank n in a charge 
(or mass) list when we arrange the clusters in the order of 
decreasing size. For instance the charge Zi of the second 
largest cluster with rank n — 2 is one-half of the charge 
Z\ of the largest cluster, the charge Z3 of the third largest 
cluster with rank n = 3 is one-third of the charge Z\ of 
the largest cluster, and so on. In the simulations of this 
work, we averaged the charges for each rank in charge lists 
of the events: we averaged the charges for the largest clus- 
ters in each event, averaged them for the second largest 
clusters, averaged them for the third largest clusters, and 
so on. From the averaged charges, we examined the rela- 
tion between the charges Z n and their ranks n. Figure [17] 
shows such relations of Z n and n for Xe with different tem- 
peratures. The histogram is the simulated results and the 
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straight lines represent the fit with Z n oc n~ x in the range 
of 1 < n < 10, where A is the slope parameter. A is 5.77 at 
T = 3 MeV. Then we increased the temperature and ex- 
amined the same relation and obtained A = 3.65 and 1.53 
at T = 4 and 5 MeV, respectively. Up to T = 5.5 MeV, 
A = 1.00, i.e., at this temperature the relation is satisfied 
to the Zipf's law: Z n oc n _1 . When the temperature in- 
creases, A decreases; for instance, A = 0.80 at T — 6 MeV 
and A = 0.56 at T = 7. The temperature at which Zipf's 
law emerges is consistent with the phase transition tem- 
perature obtained in Fig. [TH illustrating that the Zipf's 
law is also an additional signal to determine the location of 
a phase transition. From a statistical point of view, Zipf's 
law could also be related to a critical phenomenon [2129] . 
The upper panel of Fig. [TS] summarizes the parameter A 
as a function of temperature. 



10 




4 S 6 

T {MeV) 

Fig. 18. Slope parameter A of Z n to n (top) and x 2 test for 
Zipf's law (bottom) as a function of temperature for 129 Xe at 
Pf ~ 0.38po- The arrow represents the estimated temperature 
of the phase transition. Figure is taken from [65] , 



In order to further illustrate that Zipf's law is most 
probably fulfilled in phase transition points, we directly 
reproduce the histograms with Zipf's law: Z n — c/n. In 
this case, c is the only parameter, but what we are inter- 
ested in is to check the hypothesis of Zipf's law through 
a x 2 test. The bottom panel of Fig. [TBI shows the % 2 /ndf 
for the Z n - n relations at different T. As expected, the 
minimum x 2 /ndf is observed around the phase transition 
temperature, which further indicates that Zipf's law of the 
fragment distribution occurs around the liquid-gas phase 
transition point. 
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Fig. 19. Zipf plots in six different excitation energy bins for 
the QP formed in 40 Ar + 58 Ni. The dots are data and the lines 
are Zipf-law fits. The statistical error is smaller than the size 
of the symbols. 



7 Experimental Evidences of Nuclear Zipf 
Law 

7.1 NIMROD results 



In Sec. 4.2, we gave some information on critical behaviors 
for the Texas A&M NIRMOD data based on the moment 
analysis technique . Different signals of critical behavior 
coherently pointing to the same excitation energy interval 
have been shown. In this section, we will further show the 
significance of the 5-6 MeV region in NIMROD data us- 
ing a Zipf's law analysis. In Fig. [19] we present Zipf plots 
for rank ordered average Z in six different energy bins. 
The lines in the figure are fits to the power law expression 
(Z n ) oc n~ x . FigurcirUTf) shows the fitted Zipf exponent, A 
parameter, as a function of excitation energy. As shown in 
Fig. 1191 this rank ordering of the observation probability 
of fragments of a given atomic number, from the largest 
to the smallest, does indeed lead to a Zipf's power law 
parameter A = 1 in the 5-6 MeV/nucleon range. Around 
this excitation energy, the mean size of the second largest 
fragment is 1/2 of that of the largest fragment; that of 
the third largest fragment is 1/3 of the largest one, etc. 
This is a special kind of size topology of fragment distribu- 
tions, which is very different from the equal-size fragment 
distribution expected if fragments are formed through a 
spinodal instability inside the phase coexistence region 
[22I79I80I81I82I83| . This shows the relevance of using Zipf- 
plots to explore the fragment size topology. 



Y.G. Ma: Moment Analysis and Zipf Law 



13 




i - 



1 3 S 4 6 7 S 

n 

Fig. 20. Zipf law fit to the dependences of the mean charge of 
the fragment on its rank. The different symbols represent the 
multifragmentation data of different beams with an emulsion 
target. Circles and solid line represent Pb beam at 158 AGeV, 
squares and dashed line represent Au beam at 10.6 AGeV, 
star and dotted line represent Au beam at 0.64 AGeV. Data 
are taken from Ref. [84] . 

7.2 CERN Emulsion Experiment 

The nuclear Zipf-type plot has been also applied in the 
analysis of CERN emulsion or Plastic data of Pb + Pb 
or Plastic at 158 AGeV following Ma's proposal on Zipf 
law, and it was found that the nuclear Zipf law is satis- 
fied in coincidence with other proposed signals of phase 
transition [84185) . 

Dabrowska et al. have extended these studies to the 
multifragmentation of lead projectiles at an energy of 158 
AGeV [84]. The analyzed data were obtained from the 
CERN EMU13 experiment in which emulsion chambers, 
composed of nuclear target foils and thin emulsion plates 
interleaved with spacers, allow for precise measurements 
of emission angles and charges of all projectile fragments 
emitted from Pb- Nucleus interactions. The results on frag- 
ment multiplicities, charge distributions and angular cor- 
relations are analyzed for multifragmentation of the Pb 
projectile after an interaction with heavy (Pb) and light 
(Plastic - C^H^Oi) targets. A detailed description of the 
emulsion experiment can be found in Ref. [84j . 

Figure [20] shows the Zipf-type plot for charged frag- 
ments heavier than helium emitted in multifragmentation 
events of Au or Pb projectile at different beam energies. 
The values of A exponents from fits (Z n ) ~ n~ x are 0.92 
± 0.03, 0.90 ± 0.02 and 0.96 ± 0.04 for beam energies of 
158, 10.6 and 0.64 AGeV, respectively. Within the statis- 
tical errors, the values of the A coefficient are the same in 
the studied energy interval ( < 1-158) AGeV and do not 
differ significantly from unity [84j . 

Dabrowska et al. also studied the dependence of the 
power law exponent A on the control parameter m, the 



normalized multiplicity with respect to the total charge of 
spectator particles [85] . In Fig. 1217 a) are shown the mean 
multiplicity (Nf) of fragments with Z > 3 and the mean 
number {Nimf) of the intermediate fragments. The lat- 
ter are usually defined as fragments with 3 < Z < 30. In 
Fig. [21Tb) the dependence of the exponent r of the power 
fits to the charge distribution of fragments, performed at 
different ranges of m, is also given. In this analysis, the fits 
are restricted to fragment charges smaller than Z = 16. 
At small values of m, a system has few light fragments 
and the power law is steep; at large values of m there 
are basically only many light fragments leading again to 
a steep power law. At the moderate excitation energies 
where heavier fragments appear and where we expect the 
phase transition, the exponent r has its lowest value. As 
it can be seen from Fig. I21f b). the minimum r occurs 
for m values between 0.35 and 0.55. In Fig. [2TT c) the de- 
pendence of A obtained from the fits (Z n ) ~ n~ A , as a 
function of m is depicted. The exponent A decreases with 
increasing m. Between m « 0.3 and m w 0.5 the value 
of A is close to unity and Zipf's law is satisfied. This sug- 
gests that at this value of m the liquid-gas phase transition 
might occur. It has been checked that A ~ 1 occurs in the 
same region of m irrespectively of the mass of the target 
[85) . This means that the liquid-gas phase transition oc- 
curs when a given amount of energy is deposited into the 
nucleus and does not depend on the mass of the target. 
As expected, in the case of a liquid-gas phase transition, 
the previously shown maxima in frequency distributions 
of multiply charged fragments (Fig. l2"1Ta) ) as well as a 
minimum of the power law parameter r (Fig. I21f b) ) , all 
occur at the same values of m, where Zipf's law emerges. 

7.3 Some Remarks on Zipf Law 

Campi et al. pointed out that for an infinite system, Zipf's 
law is a mathematical consequence of a power law cluster 
size distribution with exponent r ~ 2 [86] . More precisely, 
both Zipf law exponent A and Fisher scaling power-law ex- 
ponent r are connected through the formula A = 1/(t — 1) 
in an infinite system assuming that the cluster size distri- 
bution is a power-law distribution. They argued that such 
distributions appear at the critical point with r ~ 2 of 
many theories, eg. various theories of cluster formation 
but also in the super-critical region of the lattice-gas and 
realistic Lennard- Jones fluids [87] . However, the experi- 
mental fragment size distribution is mostly neither power 
law distribution nor exponential distribution except for 
some special situations. Also, the nuclear system is al- 
ways a finite system, which means that the relationship 
between A and r mentioned above is not strictly valid. 
To account for finite size effects, Bauer et al. [S5] have 
evaluated the fragment probabilities as a function of their 
rank at the critical point for a finite system with fragment 
distributions obeying to a finite size scaling ansatz. From 
this analytical evaluation, where however the assumption 
is made that all fragments including the largest are much 
smaller than the source, they suggest to extend the sim- 
ple Zipf's law to a more general Zipf-Mandelbrot distri- 
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Fig. 21. (a) Mean number, (Nf), of fragments (squares) and 
mean number, (Nimf), of intermediate mass fragments (cir- 
cles) as a function of the normalized multiplicity m. Error bars 
are smaller than the size of the squares and circles, (b) Power 
law exponent, r, of the charge distribution of fragments in dif- 
ferent intervals of m. (c) Power law exponent, A, in the Zipf's 
law (see text) in different intervals of m. Error bars are smaller 
then data points. The data is taken from Ref. [85] . 



bution |89l90j . (A r ) = c(r + fc)~ A , where the offset k is 
an additional constant that one has to introduce, and A is 
asymptotically approximated as a function of the critical 
exponent r, A = l/(r — 1) of the infinite system. 

In any case, the Zipf-type plot is a direct observable 
allowing to characterize the fragment hierarchy in nuclear 
disassembly, and as such it is a useful signal of phase tran- 
sition or critical behavior. 



8 Summary and Outlook 

In summary, the moment analysis method has been in- 
troduced and some applications to nuclear multifragmen- 
tation have been presented. Since we are dealing with a 
finite nucleus rather than infinite nuclear matter, finite 
size effects must always be discussed in the model calcu- 
lations and data analysis. Experimentally, the critical be- 
havior of nuclear disassembly can be investigated with the 
help of moment analysis. The occurrence of a fluctuation 
peak which can be extracted from the moment analysis 
method can be interpreted as a signal of critical behavior. 
Using the same analysis method as for the percolation 
model, the liquid-gas universality class exponents are ap- 
proximately obtained in nuclear multifragmentation, such 
as in EOS data and NIMROD data. This would point to 
the observation of the liquid-gas critical point or second 
order phase transition. However, when we think about the 
system size dependence of critical exponents and we con- 
sider some results using Lattice Gas Model simulations 



and other related different analysis methods, it appears 
that some open questions still remain concerning the or- 
der of the phase transition. For instance, EOS Collabora- 
tion claimed that there are continuous phase transitions 
for heavier systems, namely Au and La and first order 
phase transitions for lighter system, namely Kr. On the 
other hand, NIMROD data show critical behavior, corre- 
sponding to a continuous phase transition, for the light 
system Quasi-Ar. Different conclusions are then reached 
for similar light systems. Recent systematic analyses of 
caloric curves [6|14|20] and configurational energy fluctu- 
ations [23], indicate that heavier systems may undergo a 
first order phase transition while lighter systems can prob- 
ably sustain a higher temperature, possibly even above 
the critical point, which would make the first order phase 
transition observed in heavy nuclei become a cross-over in 
lighter systems. Concerning configurational energy fluctu- 
ations, a well pronounced peak at an excitation energy 
around 5 MeV was shown in Multics, Indra, Isis and Nim- 
rod data |23j . However, this fluctuation appears monoton- 
ically decreasing in EOS data [33] ■ Thus it deserves further 
investigations. 

Scaled factorial moments and intermittency have also 
been reviewed and some examples given to show the ap- 
parent intermittency in nuclear fragmentation. However, 
some complex ingredients in experimental measurements, 
such as mixtures of event multiplicities or temperature 
fluctuations in the data can induce spurious intermittency- 
like behavior which implies that the apparent "intermit- 
tency" can not be taken as a unique signal of the critical 
behavior. Without 47r detector upgrades allowing better 
data sorting, it remains difficult to take apparent "inter- 
mittency" behaviors as a signal of critical behavior in nu- 
clear multifragmentation. 

Finally, nuclear Zipf-type plots are introduced and 
Zipf's law is proposed to be related to a phase transi- 
tion or a critical behavior of nuclei. Around the transi- 
tion point, the cluster mass (charge) shows inversely to 
its rank, i.e. Zipf's law appears. Even though the crite- 
rion is phenomenological, it is a simple and practicable 
tool to characterize the fragment hierarchy in nuclear dis- 
assembly. The 47T multifragmentation data of heavy ion 
collision at Texas A&M University and the CERN emul- 
sion/plastic data exhibit the Zipf law around the same ex- 
citation energy deposit. The satisfaction of the Zipf law for 
the cluster distributions illustrates that the clusters obey 
at this point a particular rank ordering distribution very 
different from the equal-size fragment distribution which 
may occur due to spinodal instability inside the liquid-gas 
coexistence region. To conclude, we should mention that 
all these transition signals, such as the fluctuation peak, 
critical exponents, Fisher scaling as well as Zipf's law etc 
may not be very robust individually since we are facing 
a transient finite charged nuclear system. A unique sig- 
nal can not give any definite information as to whether 
the system is in a critical point or is undergoing a phase 
transition. Only many coherent signals emerging together 
can corroborate the observation of a phase transition or a 
critical behavior in finite nuclei. 
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